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The quasi-one-dimensional S = 1 Heisenberg antiferromagnet with a biquadratic term is investi- 
gated at zero temperature by quantum Monte Carlo simulation. As the magnitude of the inter-chain 
coupling is increased, the system undergoes a phase transition from a spontaneously dimerized phase 
to a Neel ordered or spin nematic phase. The numerical results suggest the possibility of an un- 
conventional second-order transition in which the symmetry group characterizing one phase is not 
a subgroup of the other. 

PACS numbers: 75.40.Mg, 75.10.Jm 



In the search for novel quantum phases and transitions, 
various spin models with interactions beyond the ordi- 
nary bilinear couplings have been investigated. Among 
such models, the 5 = 1 Heisenberg model with bi- 
quadratic interaction is known to possess various inter- 
esting non-magnetic ground states, such as a dimerized 
and a Haldane phase in one dimension and a quadrupo- 
lar (or spin nematic) phase in two or more dimensions. 
In most magnetic systems the biquadratic term corre- 
sponds to a higher order exchange coupling and is much 
smaller in magnitude than the bilinear interactions, bear- 
ing no large effects upon the thermodynamic behavior. 
The investigation of this model at large values of the 
biquadratic term was thus mostly of academic interest. 
Recently, however, experiments 1] on optically trapped 
ultra-cold atoms opened new possibilities of flexible and 
controllable quantum systems. In particular, cold atoms 
in a periodic potential (e.g. an optical lattice), exhibit a 
superfluid to Mott-insulating phase transition 0], which 
was subsequently studied by quantum Monte Carlo sim- 
ulation [3, |j, H, |6j showing a potential capability of real- 
izing various strongly-correlated quantum systems. More 
closely related to the present study, it was shown 0, 
that Na atoms trapped in an optically generated periodic 
potential, forming a Mott insulator with one atom per 
lattice site, the hyper-fine 5 = 1 spins in the insulating 
phase may be described as an 5* = 1 Heisenberg antifer- 
romagnet with a large biquadratic term, comparable in 
magnitude to the bilinear term. 

In this Letter we present the results of quantum Monte 
Carlo simulations of the two-dimensional model with 
strong spatial anisotropy. When the anisotropy is varied 
this model exhibits a quantum phase transition from a 
dimerized phase to a spin nematic phase or a Neel phase, 
as shown in Fig. ^ The dimerized phase present in 
one dimension does not survive in the isotropic two- and 
three-dimensional model as conjectured by reference 0, 
but is stable only in anisotropic quasi-one-dimensional 
systems. Of particular interest is the quantum phase 
transition between the dimerized and the spin nematic 
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FIG. 1: A schematic phase diagram of the S = 1 Heisenberg 
antiferromagnet with a biquadratic term. Two arrows indi- 
cate values of 8 (see text for the definition) at which Monte 
Carlo simulations are carried out. 



phase, at which we find that the correlation length is 
either infinite or at least very large (£ S> 100). If the 
transition is of the second order, it must fall into an un- 
conventional class of critical phenomena 0, IToL fTT| , be- 
cause the symmetry group characterizing one phase is 
not a sub-group of the other. 

The model we consider is described by the following 
5 = 1 spin Hamiltonian: 
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where the couplings are parameterized by an angle pa- 
rameter such that J, 



(L) 



Jij cos 9, and J 
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Jij sin ( 



We focus on the region 9 < and will specifically con- 
sider two cases 9/ir = —1/2 and —0.55. While the model 
generally has the obvious SU(2) symmetry, it possesses 
an enlarged SU(3) symmetry in the case of 9 = —tt/2. 

A rough sketch of the phase diagram was presented 
previously 0] and is reproduced in Fig. ^ For the one- 
dimensional system, there are several special points at 
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which exact results are available [T3L llil [Til Il6| , and the 

complete phase diagram has been established, except for 
the possibility of an intermediate phase between the fer- 
romagnetic phase and the dimerized phase [tH ITsl , 
whose presence or absence does not have any direct con- 
sequence to the present letter though |2f)j | . In two dimen- 
sions the mean-field theory [2ll predicted the zero- 
temperature phase diagram as shown in the top part of 
Fig. which was confirmed by Quantum Monte Carlo 
simulations (2^, ■ The main focus of the present Let- 
ter is the phase transition that occurs as we go from one 
dimension to two dimensions by increasing the coupling 
between chains. We introduce the anisotropy parameter 
A such that Jy = J when the bond (ij) is along the 
chain (x) to one direction and Jy = A J if the bond (ij) 
is perpendicular to the chains (the y direction). 

Our Monte Carlo simulations are based on the loop 
algorithm with non-binary loops |25|. The system sizes 
we explored range from L x = 4 to L x — 192 with a 
fixed aspect ratio L x /L y = 4. For each data point, 
we typically perform a simulation of more than 800, 000 
Monte Carlo sweeps. The inverse temperature ranges 
from /3J — L x /2 to 2L X for the system size L x < 128. 
We have checked that the estimated values of quanti- 
ties discussed here do not depend on the temperature 
once the inverse temperature exceeds L x /J. Therefore, 
we consider that we can regard the estimates in the 
present letter as zero-temperature values. We have mea- 
sured two observables: the quadrupolar moment and 
the dimerization parameter. While the former char- 
acterizes the breaking of the spin rotational symme- 
try, the latter characterizes the breaking of the trans- 
lational symmetry. We define the quadrupolar moment 
as Q = (S x ) 2 — (S v ) 2 . As dimerization order parame- 
ter, we use D(r) = (— l) ric (Hd(r + e x ) — Hd(r)), where 
Hd(r) = diag (H rtr+ex ) , and e x is the lattice unit vector 
in the x-direction. 

We first look at the special point 9 = —tt/2, which has 
an enlarged SU(3) symmetry. In this case, the quadrupo- 
lar moment can be transformed to the staggered magne- 
tization by a 'staggered' SU(3) rotation. Because of the 
symmetry change exactly at 9 — —n/2 and our previ- 
ous numerical evidence in two dimensions, we conjecture 
that the phase boundary between the quadrupolar re- 
gion and the Neel region is vertical and straight exactly 
at 9 = — it/ 2 as depicted in Fig. ^and ends at the critical 
point (marked as 'X' in Fig. [[). Along this line, both 
the quadrupolar moment and the staggered magnetiza- 
tion are finite. In order to locate the transition point, 
one usually computes the structure factors or, more con- 
veniently, the Binder parameter (the dimensionless ratio 
of the moments of the order-parameter distribution func- 
tion). However, they are integrated quantities and, in 
the present case, are dominated by short-range behavior 
of the correlation function. Therefore, we have used as 
our probe the correlation function between two distant 
points. In particular, the ratio of the two-point correla- 
tion function of the largest distance (L x /2) to that of a 
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FIG. 2: The ratio of the two-point correlation function of the 
largest distance (= L x /2) to that of a half of it (= L x /i) for 
the quadrupolar moment (the upper panel) and the dimeriza- 
tion (the lower panel) at 6 — —iv/2. The insets are finite-size- 
scaling plots with A c = 0.1286, j/q = 2.9 (the upper panel) 
and Ac = 0.1296, j/d = 2.54 (the lower panel). The vertical 
dashed lines in the insets indicate the position of zero (the 
transition point). 

half that distance (L x /A) j2|J is useful. For the quadrupo- 
lar moment, we compute Rq(L x ) = Gq (-f-) / Gq (%■) 
where Gq(1) = (Q(le x )Q(0)) . The correlation ratios for 
the dimerization D are defined similarly. The correla- 
tion ratio for the quadrupolar moment is shown in the up- 
per panel of Fig. As expected, the curves for different 
system sizes share a common crossing point, and it is near 
A ~ 0.129. A similar analysis is applied to the dimeriza- 
tion correlation function. (See the lower panel of Fig. 

Again all the curves cross each other in the vicinity 
of a single point. The location of this common crossing 
point is A ~ 0.129 and Rd ~ 0.75. For a more system- 
atic analysis, we also apply finite-size-scaling as shown in 
the insets of Fig. [2] From the quadrupolar moment, we 
obtain A c = 0.1287(8), y Q = 1/vq = 2.9(2). The criti- 
cal exponent, t]q, can be estimated from the size depen- 
dence of the correlation function Gq(L/2) by assuming 
G Q (L/2) oc L-( d - 2 +to) (D = d + z) at the critical point 
estimated above. The results is Z + rjQ = 1.2(1). It is 
tempting to compare this result with the prediction of the 
deconfmement critical-phenomena (DCP). 0,0,0], Al- 
though it is not totally clear whether such a comparison 
is appropriate in spite of the difference in the symmetry 
of the systems (the four-fold rotational symmetry in [ij 
and [T(il ] and the two- fold symmetry in the present case), 
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the present result for 7/q is certainly larger than the value 
of the universality class of the three-dimensional Heisen- 
berg model, which is consistent with the prediction of 
the DCP. A similar analysis for the dimerization yields 
A c = 0.1288(8), yv = 2.5(2). As for r? D , we do not quote 
an estimate here because a power law behavior has not 
been observed for a large enough range of the system size 
to establish a reliable estimate. 

It should be kept in mind here, however, that these 
analyses are based on the assumption of a conventional 
critical point, which may not be justified in the present 
case. In fact, there are reasons to expect poor qual- 
ity of the finite-size-scaling plot, if the present criti- 
cal phenomenon is similar to the one discussed in Refs. 
ul EH [HI • Firstl y, t here must be two length scales (such 
as £ and £vbs in UM) diverging with different exponents. 
Therefore a good scaling plot is possible only for sys- 
tem sizes greater than the shorter (but still divergent) 
scale. Secondly, due to the existence of the secondary 
fix ed p oint (analogous to the U(l) spin liquid fixed point 
in [lj) for the renormalization flow to the VBS fixed 
point, the apparent value of the scaling exponent y in 
one phase can differ from that in the other, which could 
be fatal for an ordinary finite-size-scaling analysis. 

We next look at the case 8 = — 0.55-7T and examine the 
nature of the transition at the point 'Y' in Fig. 2] In the 
previous case of 6 = — n/2, if we regard the ground state 
as the Neel state rather than the spin nematic state, the 
symmetry group of the magnetic phase is a sub-group of 
that of the dimerized state. Therefore, the transition is 
not 'unconventional' even if it is of the second order. This 
is in contrast to the case studied previously |TotlTl| | where 
the lattice rotation (by 90 degrees) is spontaneously bro- 
ken in dimerized phase whereas it is not broken in the 
Neel state. At 9 = — 0.557T, the SU(3) invariance does not 
exist and the model has only the natural SU(2) invari- 
ance. In the spin nematic phase, the symmetry related to 
the spin space is reduced to Zi due to the spontaneous 
quadrupolar moment whereas the spatial symmetry is 
the group of all the lattice translations. In the dimerized 
phase, on the other hand, the full SU(2) symmetry exists 
in the spin space whereas the spatial symmetry is reduced 
to the sub-group that does not contain odd translations 
in the x-direction. Therefore, in the conventional picture 
of second-order transitions, the transition between two 
such phases must be of the first order. 

Recently, however, the possibility of second-order tran- 
sitions was discussed extensively and a new paradigm, 
namely the deconfinement critical phenomenon, was pro- 
posed. 0, 0, [n|, In this scenario, the VBS state is 
viewed as a state in which S — 1/2 quasi-particles are 
confined by the potential that is sustained by the spon- 
taneous VBS order to form 5 = 1 elementary excitations. 
As the quantum phase transition is approached, the VBS 
order is weakened and so is the confinement potential. 
Eventually, the quasi-particles are deconfined at the crit- 
ical point, giving rise to quasi-long-range order. It was 
also predicted that the two phases with apparently unre- 
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FIG. 3: Similar to Fig. Hbut for 9 = -0.55tt. The parame- 
ters used for the insets are A c = 0.0788, j/q = 2.9 (the upper 
panel) and A c = 0.0800, j/d = 2.8 (the lower panel). 



lated symmetries may have a direct second-order transi- 
tion due to this mechanism. For this hypothetical class 
of phase transitions, there are efforts to find an exam- 
ple that exhibits such a transition p^ . However, so far 
no phase transition of any model is confirmed beyond 
reasonable doubt to match the theoretical description of 
DCP. In particular, it was shown in Ref. j2^| that the 
DCP may not generically exist in models with U(l) sym- 
metry. 

In Fig. |3| we plot the correlation ratio at 6 = — 0.557T, 
We can see that both sets of curves cross near A = 0.078, 
an indication of a single phase transition between two 
phases. For the quadrupolar moment, the finite-size- 
scaling analysis yields A c = 0.079(1), j/q = 2.8(2). 
From the size dependence of the correlation, we obtain 
z + VQ — 1-25(20). For the dimerization, on the other 
hand, we obtain A c = 0.080(1), y D = 2.8(2). The system 
size dependence of the correlation ratio is not in perfect 
accordance with a first-order transition: if the transi- 
tion is of the first order, the correlation ratios for both 
the quantities must be 1 simultaneously at the transi- 
tion point. However, the ratio at the transition is almost 
constant and shows no visible trend towards unity. The 
correlation function decreases algebraically at the transi- 
tion at least up to a distance of 128 lattice spacings. 

Finally, in Fig. we plot the energy per site as a 
function of A. The system size dependence of the energy 
per site is small and no significant structure visible in 
this figure, in particular near the transition point A ~ 
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FIG. 4: The ground-state energy per site for 9 — — 0.55-7T, 
after subtraction of the dominant linear contribution 1.4A, 
highlighting the absence of any cusp in the energy near the 
transition point (A ~ 0.079). The inset shows the raw data 
before subtraction. 

0.079. It is still possible that the present transition is a 
very weak first-order transition of which the jump in the 
first derivative is impossible to observe within the present 
computational accuracy. Another open possibility is that 
the present phase transition is actually two second-order 
transitions very close to each other. While we cannot 
completely rule out these possibilities, we have obtained 
no evidence that positively supports either possibility. 



The last and the most interesting scenario consistent 
with the present numerical evidences is a single second- 
order transition. If this is the case, the transition is un- 
conventional as discussed above, and may correspond to 
the predicted DCP scenario which is absent in generic 
U(l) symmetric models. [3^ Because of the difference in 
the spatial symmetry of the model, namely, the lack of 
the lattice rotational symmetry in the present model, 
it is not clear how the phase transition studied in this 
Letter is related to the specific mechanism of deconfincd 
criticality that involves monopole-type excitations. How- 
even the generic scenario proposed in the previous stud- 
ies 0, [HJ, HU , in which a secondary fixed point exists and 
the translational-invariance-breaking field becomes rele- 
vant at this point, may still apply to the present case. 
While the computations shown in this Letter are already 
fairly extensive (roughly 20 cpu-years), even larger-scale 
computations are being prepared for in order to single 
out the correct scenario for the transition in the SU(2)- 
symmetric model and and to clarify its mechanisms. 
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